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Abstract 

Reliable communication over delay-constrained block-fading channels with discrete inputs and mis- 
matched (imperfect) channel state information at the transmitter (CSIT) is studied. The CSIT mismatch 
is modeled as Gaussian random variables, whose variances decay as a power of the signal-to-noise ratio 
(SNR). A special focus is placed on the large-SNR decay of the outage probability when power control 
with long-term power constraints is used. Without explicitly characterizing the corresponding power 
allocation algorithms, we derive the outage exponent as a function of the system parameters, including 
the CSIT noise variance exponent and the exponent of the peak power constraint. It is shown that CSIT, 
even if noisy, is always beneficial and leads to important gains in terms of exponents. It is also shown 
that when multidimensional rotations or precoders are used at the transmitter, further exponent gains 
can be attained, but at the expense of larger decoding complexity. 
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I. Introduction 

Temporal power control across fading states can lead to dramatic improvement in the outage 
performance of block-fading channels [1]. The intuition behind this phenomenon is that power 
saved in particularly bad channel conditions can be used in better channel realizations. Power 
control over block-fading channels was originally studied under the idealistic assumptions of 
perfect channel state information (CSI) at the transmitter (CSIT) and Gaussian signal constella- 
tions [1]. Acquiring perfect CSIT is however a challenging task due to the temporal variation of 
wireless media, as well as due to the processing and transmission delay. This motivates a large 
body of works studying fading channels under less optimistic assumptions about the CSIT; see 
for example [2], [3] and references therein. 

This work considers a block-fading channel with discrete input, where the transmitter has 
access to a noisy version of the CSI. Similarly to [4], we model the CSIT noise as Gaussian 
random variables whose variances decay as a negative power of the signal-to-noise ratio (SNR). 
Such a noise-corrupted CSIT model is well motivated and studied in the literature; see for 
example [5]-[7]. The rate of decaying of the CSIT noise can also be related to practical 
parameters in wireless systems [8]. Unlike the constant-power variable-rate scenarios, studied 
e.g. in [4], [9], we consider a power-controlled constant-rate system. In sharp contrast to the 
assumption of using Gaussian codebooks [4], [8], [10]— [13], the current work assumes that the 
input symbols are taken from a discrete distribution such as M-QAM or PSK. 

Focusing on the high signal-to-noise ratio (SNR) regime, we establish the diversity gain of 
block-fading channels under the noisy CSIT model of interest. Note that unlike in the diversity- 
multiplexing tradeoff analysis [14] where the code rate grows with the SNR, herein we keep 
the constellation size to be 2 M at all values of the SNR and we do not let the code rate scale 
with the SNR. We show that the diversity gain of coded-modulation systems can only match 
that provided by the ideal Gaussian codebooks when the ratio between the code rate and the 
constellation size is sufficiently small. The results shed some light into the interplay in the 
high-SNR regime between the number of receive antennas, the number of fading blocks, the 
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constellation size, the code rate, as well as the SNR exponent of the CSIT noise variance and 
the peak exponent constraint. 

This paper is organized as follows. The system model is given in Section [III Section [HI] 
introduces the fundamental concepts underlying our analysis. Section [IV] presents our main 
results for the outage exponent with imperfect CSIT. Section |Vl draws our final considerations. 
The proofs of our results can be found in the appendices. 



Consider transmission over a block-fading channel with B sub-channels, where each sub- 
channel has a single transmit and m receive antennas. The mutually independent channel gains 
hi, . . . , Kb have independent and identically distributed (i.i.d.) complex Gaussian components 
with zero means and unit variances. The channel gains are constant during one fading block 
but change from one block to the other according to some ergodic and stationary Gaussian 
process. This models a typical delay-limited scenario in wireless communications, where the 
delay constraint dictated by higher-layer applications prevents the system from fully exploiting 
time diversity [1]. 

The corresponding discrete-time complex baseband input-output relation for the ith sub- 
channel can be written as 



where Y i E C mxL is the received signal matrix corresponding to block i, x { £ C L is the 
transmitted vector in block i, x T denotes the transpose of x, and W \ E £; rnxL denotes the 
complex additive white Gaussian noise whose entries are i.i.d. with zero means and unit variances. 
We denote the block length by L and the power in block i by P { . Hence, a codeword corresponds 
to BL channel uses. 

We assume perfect CSI at the receiver (CSIR), i.e., the receiver has perfect knowledge about 
all the channel gains and the powers Pj. Furthermore, we assume that the transmitter has access 



II. System Model 
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to a noisy version hi of the true channel realization hi, so that 

hi = hi + e i} i = l,...,B (2) 

where e, ; G C m is the CSIT noise vector, independent of hi, with i.i.d. Gaussian components 
with zero mean and variance of. This model of the CSIT has been well motivated in many 
different contexts, such as in scenarios with delayed feedback, noisy feedback, or in systems 
exploiting channel reciprocity [5], [6]. We further assume, as in [4], that the CSIT noise variance 
decays as a power of the SNR 

al = SNIT dc (3) 

for some d c > 0. Thus we consider a family of channels where the second-order statistic of 
the CSIT noise varies with SNR. If the CSIT for example is estimated from the reverse link 
due to reciprocity, its quality will depend on the SNR of reverse link and not the forward link. 
However, while the SNRs of the forward and reverse links are different, this difference will 
be fully captured by changing the values of d c . For convenience, we introduce the normalized 
channel gains 

hi = —hi. (4) 

Given hi then hi is complex Gaussian with mean ^hi and a scaled identity covariance matrix. 
Let ji = \\hi\\ 2 be the fading magnitude of block % and 7 = [7i---7b]. Further denote 

li = \\h~1W 2 , li = \\hi\\ 2 , 7 = [7i ■ • -7b] and 7 = [71 • • -Jb]. 

The system model and CSI assumptions are summarized in Fig. \\\ 

III. Preliminaries 

We assume transmission at a fixed-rate R using a coded modulation scheme M. C C BL of 
length BL constructed over a signal constellation X C C of size 2 M such as 2 M -PSK or QAM. 
We denote the codewords of M. by x = (xj, . . . , x^) T e C BL . We assume that the signal 
constellation X has zero mean and is normalized in energy, i.e., E[X] = and E[|X| 2 ] = 1, 
where X denotes the corresponding random variable. We denote the input distribution as Q(x). 



4 



Transmitter 



■►( x 



a 

Wi 



Yb 
-►( I ) ► 



Receiver 





: h B : 

L . . r . . J 


w B 


CSIT 




CSIR 


7 







Fig. 1. System model and CSI assumptions. 



With these assumptions, the instantaneous input-output mutual information of the channel is 
given by 

1 3 



where 



Ix(s) = E 



log 2 



(6) 



is the input-output mutual information of an additive white Gaussian noise (AWGN) channel 
with SNR s using uniformly the signal constellation X. 

The outage probability is commonly defined as in [15], [16] 

P out (R) = Pr{/( 7 ) < R}. (7) 

In this work, we are interested in the SNR exponents of the outage probability [14], [17], i.e., 

(8) 



, A log Pout (#) 

dout = hm 



SNR^oo log SNR 
We adopt the notation g(SNR) = SNR a <(=> lim SNR _ logg(SNR) 



log SNR 



It has been shown in [17], [18] that the outage exponent without CSIT is given by 

d out = md sh (R) (9) 

where 



d sh (R) = 1 + 





= B- 


' BR' 


[•(-i)J 






M 



+ 1, 



(10) 



with L^J being the largest integer that is not larger than x and \x] being the smallest integer 
that is not smaller than x, is the Singleton bound on the block-diversity of the coded modulation 
scheme M [17], [19], [20]. 

Due to the availability of a noisy version of the channel 7, the transmitter can adapt the 
transmitted powers Pi to the channel conditions. In this work, we consider power allocation 
algorithms that treat the noisy CSIT 7 as if it were perfect. We consider an average power 
constraint, such that 



E 



i=i 



= E [P(7)] < SNR 



(ID 



where we have denoted P(~f) = ^ J2f=i as me instantaneous average (or normalized total) 
power allocated given the noisy channel observation 7. Thus the SNR herein has the meaning 
of the average transmit power over infinitely many fading blocks. It is well known that power 
allocation with average power constraints yields significant gains with respect to power allocation 
with peak power constraints both in terms of exponents and absolute outage probability [1]. In 
order to give a more accurate characterization of the system behavior under practical peak-to- 
average power limitations, we also introduce a peak-to-average power constraint of the form 



P{l) < SNR d p- k 



(12) 



where d pcak is interpreted as the peak-to-average power SNR exponent. The case d pcak = 1 
represents a system whose allocated power is dominated by the peak-power constraint. Asymp- 
totically, this yields the same exponent of a system with no power control. By allowing d pcak to 
take an arbitrary value, we can model a family of systems with different behavior in the peak 
power constraint. Note that in the high-SNR regime of interest, we can for example scale the 



right hand side of (1121) by a constant without changing any conclusion. That is, any constant, 
finite ratios between the peak and the average power provides the same asymptotic behavior as 

dpeak 1 • 

The corresponding minimum-outage power allocation rule is the solution to the following 
problem 

Minimize P out (R) 



< SNR 



(13) 



subject to E[iEf=i^(7) 

sEiiW) <SNR^ 

Pi(ff)>0, i = l,...,B. 
Solving this problem even numerically is difficult in general, given our noisy CSIT model and 
the discreteness of X. To date, only in the case of perfect CSIT, the minimum outage power 
control rule is known [21], along with its asymptotic behavior. The algorithm in [21] would 
actually be used in our case by a transmitter that is ignorant of the imperfectness of the CSIT. 
Nevertheless, we can characterize the asymptotic behavior of the optimal solution in the high 
SNR regime. Following the footsteps of [14], we note that the outage exponent of the optimal 
algorithm is the same as that of a power control system that allocates power uniformly across 
the blocks, i.e, Pi(j) = P (7), Vi = 1, . . . , B. This is because we can lower- and upper-bound 
the instantaneous input-output mutual information as 

B B B 

^X>0P(7)7*) < ^f>(^(7)7*) <Y,^HBP{l)li). 

i=l i=l i=l 



(14) 



Since B is a finite constant independent of the SNR, it does not change any asymptotic behavior 
of our interest. 



IV. Asymptotic Behavior of the Outage Probability 

A. Main Results 

In this section, we study the asymptotic behavior of the outage probability. In particular, our 
main results in terms of outage SNR exponents are stated as follows. 
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Theorem 1: Consider transmission at rate R over a block-fading channel described by (OQ) 
with Rayleigh fading with mismatched CSIT modeled by © with inputs drawn from X. The 
transmitter uses power control with an average power constraint (fTTI) and a peak-to-average 
power constraint (fT2l) . Then, the outage exponents are given by 



To illustrate the above theorem, in Fig. [2] we plot the outage exponents for B = 4, m = 1 
with no CSIT (or d e = 0) and with noisy CSIT with d c = 1,2 when c? pea k > 1 + md e d s ^(R). 
As we observe from the figure, increasing d c yields a better exponent. Note that in this case, 
when the CSIT is perfect the exponent is infinitely large [21]. Observe, however, that even in 
the presence of imperfect CSIT, large gains are possible by using power control, with respect to 
the uniform power allocation case. In many practical systems we typically have d c < 1 and that 
in such scenarios d c can be related to the Doppler shift [8]. In principle, achieving d c > 1 may 
also be possible by means of power control in the feedback link [11]. Note that our main result 
in Theorem Q] (and Theorem |2]) also holds for nonzero-mean hiS (Rician fading), because the 
asymptotic diversity gain only captures the slope of the outage probability, which is the same 
for zero and nonzero-mean hj's. 

To get some insight into the problem, let us take a closer look at the results of Theorem \T\ 
in some special cases. In the extreme case <i pea k = 1, which implies that the average and 
peak power have the same exponent, we obtain d(R, d e , 1) = md s b(R), which is the outage 
exponent for a system with short-term power control [21], or no power control [17]. Since a 
system with short-term power constraints cannot allocate power across multiple codewords, it is 
logical that the resulting outage exponent is independent of the quality of CSIT. Increasing <i peak 
subsequently leads to an improvement in the outage performance. However, when <i pea k exceeds 
a certain threshold, there is no extra diversity gain by increasing a? pea k further (the diversity gain 




md sh (R) (1 + md sh (R)d e ) <i peak > 1 + md sh (R)d e . 




(15) 



Proof: See Appendix lAl 



s 




Fig. 2. Outage exponents for B — 4, m — 1 and d pca k > 1 + md c d s ^{R). 



is "saturated" due to the limitation on the accuracy of the CSIT). In other words, a stringent 
constraint on the peak power exponent leads to a lot more pronounced detrimental effect in the 
case of accurate CSIT (large d e ) than in the case of very noisy CSIT (small d e ). 

In the limiting case d e j 0, i.e., very noisy CSIT, we have d(R, d e , <i pea k) — > mdsb(R), which is 
again exactly the outage exponent when there is no CSIT [17]. In this case the outage exponent 
is also independent of rf pe ak 5 because the transmitter always uses a constant power in the order 
of SNR 1 . The case d e I also represents the scenarios in some practical systems in which the 
CSIT noise variance does not decay with the SNR. If the CSIT noise variance has such an "error 
floor" in the high-SNR regime, then no extra diversity gain can be obtained from power control. 



9 

On the other hand, in case d e — > oo, i.e. when the CSIT noise variance decays exponentially 
or faster with the SNR, then d(R, d c ) — > oo, VR < M, as long as the peak exponent constraint 
is also relaxed to satisfy d pcak > 1 + md s ^(R)d c . For strictly positive and finite d c , using power 
control, even with noisy CSIT, provides an extra diversity gain of (md sh (R)Yd e compared to 
the no-CSIT case, as long as the peak power constraint is sufficiently relaxed. The presence of 
the factor m 2 also parallels with the diversity-multiplexing tradeoff result obtained in [8] for 
MIMO channels with Gaussian inputs. 

We also learn from the analysis in Appendix [A] that at high SNR, when d pea ^ is sufficiently 
large, the dominant outage event occurs when exactly —1 of the channel gain estimates ^'s 
are much larger than the noise variance erf, and the remaining B — + 1 channel estimates 
have the same order of magnitude as of. For example, when the rate is sufficiently small such 
that BR < M then a typical outage event occurs when all B channel estimates are in the order 
of the CSIT noise variance, leading to the maximum diversity gain of mB(l + mBd c ). When 
^peak is sufficiently small, however, the system cannot "invert" the worst channel realizations and 
the peak exponent becomes the limiting factor. For example when rf pea k < d e then the dominant 
outage event happens even when all the channel estimates are very accurate (significantly above 
the CSIT noise level). 

B. Improving the Outage Exponent with Rotations 

In [22], it is shown that a simple precoding technique can be used to improve the outage 
exponent over fading channels with discrete inputs and uniform power allocation. In this section, 
we demonstrate how the idea in [22] can be applied in the current noisy CSIT setting of interest 
to further improve the outage exponents. In order to avoid cumbersome notation and to simplify 
the presentation, we remove the peak exponent constraint (setting <i pea k = oo), focusing only on 
the effects of the CSIT noise. 

In the following we briefly recall the precoding technique of [22] . First consider reformatting 



to 



the codewords x £ M. as matrices 



X 



£ C 



BxL 



(16) 



We now obtain X as 



where 



X = MS 



'Mi ^ 



(17) 



M 



£ C 



BxB 



(18) 





yo o m k j 

is a unitary block-diagonal matrix, and the entries of S £ C BxL belong to the signal constellation 
X with size 2 A/ symbols. The matrices Mi, . . . , Mk £ (£ NxN are the K unitary rotation 
matrices of dimension N each. Fig. [3] illustrates the above construction. These rotation matrices 
are required to have full diversity, i.e., 



M k {s - s') ± 



(19) 



componentwise, for all x ^ x' £ X N . This implies that if the vector (s — s') has a positive 
number of nonzero entries, then, its rotated version will have all iV entries different from zero. 
The reader is referred to [22] for more details on the construction and to [23] for a detailed 
discussion on the design of full-diversity rotation methods. 

According to [22], with no CSIT we obtain the following exponent 



d out = mdX(R) 



where 



<C(R) = N 1 + 



B 

N 



1 - 



R 

M 



B + N-N 



BR 

MiV 



(20) 



(21) 



With noisy CSIT, completely similarly to the previous section we have the following result. 
Theorem 2: Consider transmission at rate R over a block-fading channel described by (OQ) 
with Rayleigh fading with mismatched CSIT modeled by © with inputs obtained as the rotation 
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Fig. 3. Code construction with rotations. 



of a coded modulation scheme over X as described by (fTTT) . using full diversity rotations. The 
transmitter uses power control with an average power constraint (fTTj) . Then, the outage exponents 
are given by 

d(R, d e ) = md$(R) (l + md^(R)de) . (22) 



Proof: See Appendix O ■ 
We illustrate in Fig. @] the effect of full-diversity rotation matrices on the outage exponent of 
the coded modulation system with mismatched CSIT. This precoding method clearly leads to a 
higher diversity gain even at high code rates, at the expense of increasing receiver complexity. 

In the special case N = B, i.e. when a single matrix that rotates all B output symbols is 
used, then d(R,d e ) = mB(l + mBd c ). This is the maximum diversity gain we can achieve in 
this scenario, even with codes drawn from a Gaussian ensemble [8]. For a large N, however, 
the receiver complexity will increase exponentially, as this rotation will require joint decoding, 
taking the output of blocks of N sub-channels into account. Note also, that, since this strategy 
yields the optimal exponent, in terms of exponents, there is nothing to gain in optimizing the 
full precoding matrix. Using power control and a full-dimension full-diversity rotation matrix is 
sufficient. 
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Fig. 4. Outage exponents for B — 4, m = 1, ck = I and full-diversity rotations of size N = 1 (dotted line), N — 2 (dashed 
line) and N = 4 (solid line). 



V. Conclusion 

We have studied the asymptotic behavior of the outage probability for code modulation over 
block-fading channels under the assumption that the transmitter has access to a noisy version of 
the instantaneous channel gains. We showed that power control even with mismatched CSIT is 
still very beneficial in improving the outage performance of the system. Our results shed some 
light into the interplay between different parameters in a coded modulation system, including 
the constellation size, the code rate, the quality of the CSIT, and the peak power requirement. 
Determining the outage exponents in a more general multiple-input multiple-output remains an 
interesting open problem. 
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Appendix A 
Proof of Theorem Q] 

Since we are interested in the high-SNR regime, let us invoke the standard change of variables 

lo g7i anr | A. — lo S7» 
logSNR u " l logSNR - 



as in [14], a t = — . log s X' R and 6>i = — , los JL • We also perform the change of variable 7r(7) 



4 Ml) 

V / log SNR " 

The power constraint (fTT|) asymptotically becomes [8], [24] 

y SNR" ( ^/(7)^7 < SNR 1 . (23) 

Notice that the 7/s are mutually independent and follow Chi-square distribution with 2m degrees 
of freedom. Also, we have E[ 7i ] = E[||h,j|| 2 ] — E[||ej|| 2 ] = SNR . Changing variables from 7 to 
a, we readily obtain 

/ SNR^ (S) SNR~ m &i da < SNR 1 . (24) 

Herein we have neglected the terms irrelevant to the SNR exponent, noticing that for any 
set containing on < 0, its probability measure decays exponentially in SNR [14]. Applying 
Varadhan's integral lemma [25] we then have 

sup < 7r(a) — m^^dj > < 1. (25) 

Since outage probability is a non-increasing function of transmit power, we conclude that with 
the optimal power allocation, 



7T 



a) = min \ d peak , 1 + m ^ a h j (26) 



where we need to introduce <i pea k to take into account the peak constraint (U2l) . 

From [17] it is known that as SNR — > 00 the mutual information in sub-channel 2, (P(7) 7 j), 
tends to either M or depending only on the behavior of the term 

P(7) 7 j = SNR^SNR^SNR - ^ = SNR min ( dpcak ' 1+mE ?= 1 " J )" dc ~°\ (27) 

In particular, if oii < n(a) — d c then 1% (P(y)'Ji) —* M bits per channel use. Otherwise 

Ix {P(l)li) - 0. 
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Thus the asymptotic outage set is given by 
O 



ja, dc : ^ 1 < min ^d peak , 1 + m^ct,^ - 4^ < ( 28 ) 



where l(-) is the indicator function. We then have 

Pout(iZ) = / /(7l7)/(7)d7d7 



(29) 

/(a|S)/(a)(iQ;(iS. 



Notice that /(7I7) = Ylf = i f(ji\li), where the conditional p.d.f /(7i|7i) is a non-central chi- 
square one with 2m degrees of freedom. In Appendix |B] we asymptotically expand the integral 
(1291) . showing that the outage exponent is eventually given by 

d(R\ 4, (ipeak) = min(4, ...,d B ) (30) 
with d n being defined such that 

„ B-n B 

I \\ §NR- m& *- m ^ \\ SNR~ m&: > dotdet = SNRT* 1 (31) 



'OnBn i=1 j =B _ n+1 



where 



B n = {a, a : {a x > 0, a x > d c } f] ■ ■ ■ H > 0, & B - n > 4} 

n{0 < acB-n+i < 4, a B -n+i = oiB-n+i - 4} Pi • • • fl {0 < a B < d e , a B — a B — 4}} • 

(32) 

Thus applying Varadhan's integral lemma [25] gives 

( B B-n \ 

d n = inf < m > dj + m > a, > . (33) 



i=l i=l 



Recall from <|28]) that 

n / / n \ \ 



: ^ 1 < min ^ pea k, 1 + m ^ % j 



(9 ■ <; a. a : > 1 ( n, ^ : nun j r/,„, lk . 1 + /// > n :j j - d c j < ^ 

Over B n , we have that «j = &i — d e for all i > B — n + 1, thus 

O = < 6c, ac : ^ 1 I «i < min I ^peak, 1 + m ^ dj ; ) - d c 

/ / /, \ \ BR 



^2 1 I &i < mm I rf peak, 1 + m ^ 
j =J B-n+l \ \ j=l 



M 
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To compute d n , we consider two mutual exclusively cases. 

Case 1: d pcak < 1 + mj] . =1 aj. We denote the SNR exponent over the intersection of this 



region and B n as dn . Then 

B-n 



B 



O 



CK, S : ^ 1 («i < dp Cak -d e ) + ^ 1 (<*» - ^Peak) < 



BR 



(35) 



i=l i=B-n+l 

Case 1.1: If d pcak < d c then 1 («; < <i peak — d c ) — 0, Vz G {1, . . . , -B — n}. The outage set 
reduces to 

U BR 



O 



OL -. ^ 1 («i < ^pcak) < 



M 



(36) 



i=B-n+l 

Because for i = 1, . . . , B — n, the terms &i and «j are not present in the outage set, we have the 



optimal solution to (1331) a* 



a 



* B _ n = and J2f=i n «i = max(d peak - 1, m(.B - n)rf e ), 



due to the constraint <i poak < 1 + m J2f =1 dc j. 



There are only n terms in the summation in (1361) . thus if n < ^ then 

c?W = max (d peak - 1, m(.B - n)d c ) . 
But since n < ^ < B and G? peak < c? c we have m(B — n)c? c > m(B — n)d peak > d peak — 1. 



Thus 



if n < ^ 



djj) =m(B-n)d e 



If n > then without the constraint <i peak < 1 + m ^f =1 a? we readily obtain the solution 



to 03b: a 



B-n+l 



d£ = 0. Taking the 



constraint m $^,- =1 > ^poak — 1 into account we have 

d£ ] = max ^d peak - 1, m(B - n)d e + md peak | n 
But c? pcak < d c thus we have 



~M 



+ 1 



(37) 



m( J B - n)d c + md peak 
and 





'BR' 




M 



+ 1 > md 



peak 



dn } = m ( 5 - n K + "l^peak ( Tl 







(- 




M 






^72 — 






M 



+ 1 > md peak > d peak - 1, 



+ 1 • 



(38) 
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In summary, if d pea ^ < d c then 



m(B - n)d e if > n, 

<#> = < (39) 

m(S - n)(i e + md peak (n - \^f \ + l) if ^ < n. 
Case 1.2: On the other hand, if rf pe ak > rf e > then for i = B — n + 1, . . . , B we have < 
dpeak) = 1 because in B n , 6a < d e for these values of i. The outage set reduces to 

{ B ~ n BR ] 

a, a : 1 (a, < rf pcak - d e ) < — - n\ . (40) 

If if < n then d£° = oo because the set of "bad" channel realizations is empty [13]. Intuitively, 
in this case we have access to "perfect" knowledge about n channel gains which we can then 
use to successfully "invert" the channel gain (since d peak is sufficiently large and does not pose 
any restriction). Consequently we can achieve exponential decay in the outage probability for 
all rates R < Mp. 

If > n then, due to the total absence of ctj in (|40l) . the optimal solution to (1331) satisfies 
J2f=i&i = max (^peak — l,m(B — n)d e ), where we have taken into account the constraint 
<^peak < 1 + m J2f=i®j- As f° r ^«' S ' f rom <l40l ). we see that at the optimum points, there are 
exactly — n\ — 1 of the a^s that are equal to zero, and the remaining B — n— \^ — n \ +1 
variables are all equal (or arbitrarily close to from above, strictly speaking) to d pca k — d c . 

Finally we have 

( BR 1 \ 
B - n + 1 - -r— - n j + max (d pcak - 1, 771(5 - n)d c ) . (41) 

In summary, if d pea k > <4 then 

i(d peBk - d e ) (B - n + 1 - \j^-n])+ max (d poak - 1, m(B - n)d e ) if ^ > n 

BR , 

(42) 

Case 2: d pca k > 1 + rra X/f=i Note that over £> n we have Y^=x &j > (B — n)d c thus Case 
2 can only happen if d peak > 1 + m(B — n)d e . For n such that d peak < 1 + m(£> — n)d c , we use 



dl 1 ) 



my 

if^<//. 
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the convention dff = oo. Then, over B n 
O = < a. a 



B-n / B \ B / B \ 

: 1 I aj < 1 + m Q!j — <i c J + 1 ( < 1 + m cfe^ J 

i=l V i=l / i=B-n+l V i=l / 

<^ a, a: : ^ 1 I «i < 1 + m a J : _ " c I < ~ 



< 



BR 



(43) 



Again if 4^ < n then the outage event decays exponentially in the SNR. We readily obtain 



ft' -i 



a 



a: 



B 



0. We also have a* = 1 + m(_B — n)d e — d e , 



for exactly B — n — — +1 of the cn/s, and the other c^'s are zero. Thus 



m(5 - n)d e + m (B - n - \ £j£ - n] + l) (1 + m(B - n)d e - d e ) if ^ > n, 

BR 



OO 



\f on < 



We now combine the results in Case 1 and Case 2 to find the outage exponent 



d(R, d e , (ipeak) = min(d , ••-,d B ) = min(4 , d!£\ . . . , d^',d^'). 



7(1) jph 



(44) 



(45) 



If ^peak < <4 men me dm s are given by (|39l) . Furthermore, we have rf pca k < d c < l + m(B — 

r<y\ 

n)d e , Vn = 0, . . . , B — 1 thus <i„ = oo for these values of n. For n = B, we also have from 
(1441) that dn = oo. Thus in this case 

d(R, d e , rfp eak ) = min(4 1) , df\ ...,d% ] ). 

From (l39l) we have 



(46) 



Also from (1391) and from the fact that rf pca k < d c , we have 



^ f BR ~\ > ■ ■ ■ > — ^-^peak 







(- 




M 



+ 1. 



(- 


'BR' 




M 



+ 1 . 



(47) 



(48) 



Thus we finally have 

d(R, d e , (ip e ak) = min ( md e ( £ - 

™d p eak ( B 



BR 



M 



+ 1.1, md 



peak 



B 



BR 



+ 1. 



BR 



+ 1. 



(49) 
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The analysis also reveals that the dominant outage event occurs in the region Bb, i.e., when all 
the channel estimates 7i's have a much large order of magnitude than the CSIT noise variance. 
More specifically, in the typical outage event, B — \jjf~\ + 1 of the channel estimates are in 
the order of SNRT dpcak , canceling out the maximum power that can be allocated to any channel 
realization. Thus the limiting factor in this case is the peak exponent <i peak . 

We now consider the case d pcak > d e , where the di^'s are given by (|42|) . There are three 
possibilities. 

Case A: If d pcak > 1 + mBd e then <i peak > 1 + m(B — n)d c , Vn = 0, . . . , B. Thus 
i(dpeak -d e )(B-n + l-\^-n\)+ rf pcak - 1 if ^ > 



n 



if BR< n 

11 M — n - 



(50) 



But since d pcak > 1 + m(B — n)d c 

m(d peak - d c ) I B - n + 1 - 



BR 



n 



^pe a k 1 

BR 



M 



n 



+ m(B — n)d c 



> m(l + m(B - n)d c - d c ) yB - n + 1 
But the right hand side is exactly the value of dn in (1441) when ^ > n. We conclude that 

d(R, d e , d peak ) = min(4 2) , d[ 2 \ . . . , d^'). 

Furthermore, from (l44l) we have that drg R -, < ■ ■ ■ < df < d^f 1 < oo = d^ BR -, 
Hence 

d(R, d e , c? peak ) tlr gfi i 

I M I 1 

"Si? 



(2) 



T¥l-1 

m ( I + m [ B 
'BR 



+ m\^B 
ml B — 



M 

+ 1)4 





BR" 




'BR 


BR" 




M 


+ 2- 




M 


+ 1 


M 





'BR' 


(- 




M 



+ 1 )d 



AI 
BR 

~w _ 

(51) 

The dominant outage event occurs when exactly S — + 1 of the channel gains have the 
same order of magnitude as the CSIT noise variance. The peak exponent constraint is not the 
limiting factor in this case. 
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Case B: 1 + md e (B - + l) < d pcak < 1 + mBd c . This implies ^ > \^] - 1 > 



B 



^pcak 



-i 



; . For any integer n such that n < B 
for these values of n, we have 



=m{d 



peak 



<L) [ B - n + 1 



^peak 1 



BR 



then c/peak < 1 + md e (B — n). Thus 



n 



md e (B — n) 



and dffi 



oo. 



As for n 



-'■peak 



-i 



, . . . , | ^1 — 1, then rfpeak > 1 + md e (B — n). This is similar to 



Case A, i.e., we have dn > dffl- 



c/(i?, <i e , tip 



Thus in Case B 

peakj " r° ' ■ ■ • ' V%^1 

due to the fact that d?L = • • • = d# = oo for any fe. It is readily verifiable that d!^' > 

I ~M~ I 

d r d_„,.-n > d pcak -i -| > - - ■ > ^raj 



any fc. It is readily verifiable that 

d , _n > d^ d ,_,-,>•••> <ir 2 l R -| and thus 



> 



d^R., d e} C?p ea k) (jp -B-R l 1 





'BR' 






M 



1 ) (l+m[B- 



BR 



\ \ (52) 

+ 1J4J • 

gains have the same 



The dominant outage event also happens when B — \^jt~\ +1 of the channel 
order of magnitude as the CSIT noise variance. 

Case C: d pea k < I + md e (B - + l). This implies \^] - 
any integer n such that n < ^ then n < B — dp ^ -1 leading to G? pca k < l + md c (B — n). Hence 
from (|42l) we have 



KB- l h^l_ Thus for 



n 



. m(d pcak -d c ) (B-n + 1- \%r-n\) + md e {B - n) if ^ > 
= <! (53) 

BR 



OO 



if < n 
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Since n < 4^ leads to d pca k < 1 + md c (B — n), we also have dff = oo, Vn. Thus 



(i) 



4, (ip eak ) = nun(4' 

r#i-i 



5 'y 









' BR 


'BR 






+ 2- 






+ 1 




M 




M 


M 



+ md e ( /J 
md peak ( B 



BR' 
'BR 



+ 1 
+ 1 



(54) 



Again the dominant outage event occurs when exactly B — \ B 3'\ + 1 of the channel estimates 
have the same the order of magnitude as the CSIT noise variance. Unlike in Case A and Case B, 
in this case the peak exponent rf pe ak is too small and becomes the factor preventing the system 
from achieving its full potential. ■ 



Appendix B 
Asymptotic Expansion of (1291) 

In this appendix we review the asymptotic expansion of the joint p.d.f. in (|29l , a result derived 
in [8] for the case of a single fading block. In particular we would like to study the high-SNR 
behavior of 

B 

Pout(R) = / X\f{li\li)diid%. (55) 

where /(7i|7i) is a non-central chi-square p.d.f with 2m degrees of freedom and non-centrality 
parameter ^ = SNR~ Qi+de . Changing variables to a. and a gives 

Pout(R)= [ n e - SNR ~%- SNR ~ ( " ! ~ dC V SNR ^ 
J° i=l 

x SNR s r 1(fii - ai -* ) - Si - mai J m _ 1 (sNR de ""2" a ' ; \ doLidoti. 
For each i £ {1, . . . , B}, let us define the set 



(56) 



— ati : d e — &i — ati > 0} 



(57) 
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and its complement 

A { P = {da, a.i : d e - &i - a» < 0} . (58) 

Firstly, consider the region A{ \ i.e., d e — &i — 6ti > 0, for some i. Then SNR dc_Q ^ _ai — > oo 
as SNR — > oo. But for real x > we have [26, Sec. 9.7] 

/ m _ 1 (x) = -^(l + 0(l/x)) (59) 

yZTTX 

/ de-aj — &i \ d c ~a i -a i d c- a i- a i 

thus J m _i (SNR 2 j = SNR * e ;,JNK . Grouping the exponent terms inside the 

integral (1561 ) gives 

exp (-SNR-" 1 - SNR" (&i ~ 4) + SNR ^'"'' ) exp (-SNR" &1 ) . 



Note that 



max(-oti, -(dii - 4)) > — — y — — (60) 



for any on, oci with the equality occurring iff a» = dj — d e . Therefore if ctj 7^ d; — c? e then 

— SNR~ Si — SNR"^ - ^ + SNR dc ~"'~"' = _SNR max ^~ ,5 "~^ i ~' ic ^ 

But we are considering where d e — at$ — on > 0, so max(d e — dj, —ai) > 0. Thus if 
a.i 7^ &i — d e then the outage probability decays exponentially in SNR. 

If OLi = &i — d e then the condition d c — dj — dj > leads to dj < d e . We also have 
SNR" 1 = SNR" l SNR- dc or j { = Thus we can write 



/ . m 9if{li\li)f{li)dlid%= I 9if{li)dli 

JOnA)' J Or\\& i «L,ai=6u-de\ 



Or\{ai<d e ,ai=ai-d c } 

i e- SNR ~* i SNR- fn * i da < (61) 

On{ai<d e ,ai=&i-d c } 

gSm-^dai. 



1 Or\{0<a i <d e ,a i =a i -d c } 

Herein we have denoted 

B 

Qi = j [ f(aj\&j)dajd&j. (62) 
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Secondly, consider the region A\ where d e — ct{ — &i < and thus the asymptotic form of 
the modified Bessel function of the first kind J m _i(x) with x j gives 

I m _i (SNR — = SNR (m - 1} — if— 1 . (63) 

We can then constrain «j > and dj > c/ e , because otherwise the outage probability decays 
exponentially. Thus (cf. (1561) ) 



/ gif{li\li)f{li)dlidni = [ giSmr^-^doidoi. (64) 

JOnA^ J On{ai>0,ai>d c } 

Recall that gt collects all the terms that are independent of ctj and d,. 

Thus in the asymptotic expansion of the outage probability, we need to consider 2 B regions 
nf =1 Af DO where q = 0, 1. The slowest decaying terms among these 2 B regions will determine 
the outage exponent. However, due to complete symmetry, we can assume without loss of 
generality that 6>i > ■ ■ ■ > a B . Then the number of regions need considering reduces to B + 1. 
In particular for each n E {0, . . . , B} we need to find the exponent d n where 

I f{ct\a)f{dt)dctda = SNR~ dn (65) 

J Or\{&0>--->&B-.n>dc>aB-n + l>--->6lB+l} 

with the convention a = oo and a B +i = — oo. Then d(R, d e , <i pea k) = min(<io, ■ ■ ■ , d B ). 
From (loTI) and (|64|) . we have 

/ ( a. | a ) f ( a. ) da. den 

On{iio>->a s _„>4>«B_„ + i>->a fl+ i} 

„ B-n B (66) 

/ TT SNR- m&i ~ m ^ ] [ SNR- m& idada 



OnB n i=1 j=B-n+l 



where 



B n = {a, a : {a x > 0, 6c x > d e } H • • • H {a B _ n > 0, a B - n > d e } 

n{0 < aiB-n+i < d c , a B -n+i = ois-n+i - d e } n • • ■ n {0 < a B < d e , a B — a B — d e }} . 

(67) 
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Appendix C 
Proof of Theorem [2] 

Similarly to the previous proof we have, 



d(R, d e ) = min((io, • • • , d 



b 



(68) 



where 



OG 



BR 

MA" 



if \tP\ > ... . 

r/„ = < (69) 

^ miV (f - [f 1 - K n ) (1 + m{B - n)d e - d e ) + m{B - n)d c if [f ] < ||. 

Herein 





~ n ~ 




MN 


N 





- 1. 



(70) 



In this case d n is the dominant outage exponent conditioned on the event that there are exactly 
n channel gain estimates % having a larger order of magnitude than the CSIT noise variance of. 
Note that by definition di = oo, Vi : > jjjf. Intuitively, when at least i channel gains are 
known (asymptotically) noiselessly at the transmitter, then using power control we can always 
transmit MN x [~-^] > BR bits with exponentially decaying error probability. This is because at 
worst, these known channel gains belong to the least number of rotation groups, which is | . 
Then from the definition of d n we have d N s r bri _ x \ < d N r r bk -i _ x \ 1 < • • • < d and 
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d N ( I br] _i)+i = • • • = d B = oo. Thus 
4) = min(4, ...,d B ) 



d 



Him-') 



= m^B-N ^ 


BR' 

MN 




■)) 


+ miV 


' BR' 

MN 




= m ^5 + iV - iV 


' BR' 

MN 




+ miV 


' BR' 

MN 


+ 2- 


= m ^5 + iV - N 


'BR' 

MN 




fB 
+ mN i^-- 


' BR' 

MN 


*■)( 


— m + N — N 


'BR' 

MN 


)(■ 



4 



5i? 



l + m B-JV 



Bi2 

MiV 



- 1 4 - d c 



MN 



BR 

MiV 



+ 1 



1+m B+N -N 



BR 

MiV 



BR 

Jin 

BR 

MiV 



d e - d c 



d. 



(71) 



d, - d P 
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